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We model the evolution of the galactic bulge and of the bulges of a selected sample 
of external spiral galaxies, via the multiphase multizone evolution model. We address 
a few questions concerning the role of the bulges within galactic evolution schemes 
and the properties of bulge stellar populations. We provide solutions to the problems 
of chemical abundances and spectral indices, the two main observational constraints 
to bulge structure. 
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1 INTRODUCTION 

Bulges are the spheroidal component of spiral galaxies. Mod- 
els for their formation and evolution have been dominated 
in the past by the strong parallelism with elliptical galax- 
ies, since several studies on Galactic bulge stars (e.g. Whit- 
ford 1978; Frogel & Whitford 1987; Rich 1988: Frogel 1988) 
claimed that their photometric properties were analogous 
to early-type galaxies and only old and super-solar metal 
rich stellar populations seemed to exist in bulges. Recently, 
younger stars were found in Baade's window (BW, Holtz- 
man et al. 1993), suggesting a star formation history differ- 
ent from the one attributed to elliptical galaxies, i.e. strongly 
concentrated in the early evolutionary stages and turned off 
after the first Gyr. Chemical abundances in the Galactic 
bulge have been revised (McWilliam & Rich 1994; Minniti 
et al. 1995; Ibata & Gilmore 1995a, 1995b; Houdashelt 1996; 
Frogel, Tiede & Kuchinski 1999, hereinafter FTK99) and at 
present the iron relative abundance peaks at ~ —0.3 dex. 
A similar result has been obtained by Bacells & Peletier 
(1994) and dejong (1996a,b) for external spirals, by estimat- 
ing from new photometric data a sub-solar mean metallicity 
< [Fe/H] >~ —0.2, —0.3 dex, suggesting that all bulges are 
very similar in colors and metallicities. Kinematic observa- 
tions show that bulges are rotating much more rapidly than 
the bright elliptical galaxies, suggesting that bulges are sup- 
ported by rotation at variance to ellipticals. 

The main scenarios of spiral galaxies formation con- 
straint the bulge formation: in the dissipative collapse model 
by Eggen, Lynden-Bell & Sandage (ELS, 1962) the evolu- 
tion proceeds inwards and metallicity, age and kinematics 
are correlated. Zinn (1985) proposed a merger scenario, with 
a halo formed by accretion of small fragments, justifying 



the absence of a metallicity gradient in the halo. The cor- 
relation of bulge and disc scale lengths and the similarity 
of inner disc and outer bulge metallicities stimulated alter- 
native models, such as a bar producing a concentration of 
stars in the center (Pfenninger & Norman 1990; Pfenninger 
& Friedli 1991) or a nuclear starburst ejecting a giant gas 
bubble from which stars form (Wada, Habe & Sofue, 1995). 
Wyse (1999) demonstrates that instabilities of purely stellar 
discs cannot form bulges and that the metallicity distribu- 
tion in MWG is not consistent with a bulge built up from 
accretion of satellite galaxies or clusters. 

The ELS scenario seems to recover, supported also by 
recent observations (Minniti, 1996a,b), and it suggests that 
the bulge formed through a dissipative collapse, resulting 
younger than the halo; FTK99 with IR data show that 
the metallicity presents a strong radial gradient in the in- 
ner bulge (< 1.5 kpc), in agreement with previous studies 
(Terndrup, Frogel & Whitford, 1990; 1991 and Minniti et 
al., 1995). Zinn (1993) found that the inner old halo clusters 
(R < 6 kpc) show the signature of the dissipation, and that 
there is a smooth transition from the halo to the disc glob- 
ular clusters, suggesting that they proceed from the same 
collapse. Within this revival of the dissipative collapse, Van 
den Bosch (1998, 1999) and Elmegreen (1999) discussed the 
bulge formation as due to baryonic material from halo or 
protogalaxy settling in the center of the spheroid in a short 
time. The natural trend is to form a bulge because the an- 
gular momentum is low in the center, and it results older 
than the stellar disc. In particular, Elmegreen evaluates the 
threshold density for the star formation in turbulent condi- 
tions due to the gas accretion to form the bulge: a starburst 
phase is established in a short time scale. In both cases, the 
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formation and evolution of bulges are different from ellipti- 
cal galaxies. 

In the present paper we wish to address the following 
questions; a few have already challenged the modelers and 
others have not been previously discussed in the literature: 

• Is it needed to require the merging and/or accretion of 
external material to reproduce the main observed character- 
istics of bulges? 

• Is it possible to reproduce the observed properties of 
bulge stellar populations? 

• How far goes the analogy with elliptical galaxies? 

• Which are the ages of bulges within the whole galactic 
evolution? 

To afford these questions, we adopt the multiphase evo- 
lution model, which assumes a dissipative collapse of the gas 
from a protogalaxy or halo to form the bulge and the disc, 
in ELS scenario. The model has been applied successfully to 
the Galaxy and external galaxies (see papers by Ferrini and 
coworkers) . 

Molla & Ferrini 1995 (hereinafter MF) applied first the 
multiphase model to study the evolution of the Galactic 
Bulge, considering it as a natural extension of the disc, any- 
way with a separate evolution: supernovae winds and no ra- 
dial flows from the disc were assumed. At variance with other 
models present in the literature, MF succeeded in reproduc- 
ing the metallicity data, revised at that date by McWilliam 
& Rich (1994). In the present approach, first we modify the 
modeling of the bulge, by considering a core population in 
the central region. In this scenario, the nuclear population 
is distinguishable from the bulge one as far as metallicity 
and dynamics are concerned: nuclear starts are more metal- 
rich and have higher a elements abundances. Furthermore, 
their motion must be supported by rotation, as predicted 
by dissipative collapse models. Indications about a different 
kinematics of stars in the galactic nucleus come from rela- 
tively large data sets of OH/IR stars (Lindquist et al. 1992) 
and ordinary M giants (Blum et al. 1995) in the 5-200 pc 
region. These two samples appear to trace different popu- 
lations: while the M giants appear to joint smoothly with 
the kinematics of other bulge tracers at larger radii, OH/IR 
stars have a lower velocity dispersion (oonjiR — 50 — 100 
km/s against <jm > 100 km/s) and a very high rotational 
velocity (V ro t > 100 km/s at a distance of 100 pc from the 
galactic center), against an almost negligible rotation of the 
M giants. The OH/IR stars may be associated with the nu- 
clear bulge star population. 

We extend the discussion to a set of external galax- 
ies, whose discs have been extensively studied by our group 
(Molla, Ferrini & Diaz 1996, 1997). 

The theoretical model is presented in §2; in §3 are the 
result of the new bulge model. In §4 we analyze the chemical 
evolution results for all the bulge sample, by including the 
study of possible trends with the galactic properties (Hubble 
type, total mass). In the section §5, we predict the values 
for the Lick spectral indices Mg2 and Fe52 and compare 
these predictions with observational data. Discussion and 
conclusions are in §6. 



2 THE EVOLUTION MODEL 

In the multiphase approach the star formation is considered 
a two step process: the diffuse gas forms molecular clouds, 
then stars form from cloud-cloud collisions (spontaneous star 
formation) or by the interaction of massive stars with molec- 
ular gas which induces star formation. We discussed exten- 
sively this approach in previous papers (e.g. Ferrini et al. 
1992) and fixed the relevant parameters which define the 
coupling between the phases and the zones. Interesting to 
our present paper is essentially the collapse time scale t co u, 
which defines the infall gas rate from the halo to the sec- 
ondary disc region or bulge. We already discussed how this 
varies with the galactic mass and with the galactic radius 
(Ferrini et al. 1994). Spontaneous star formation and cloud 
formation rates depend on the radial distance, too, and on 
the Hubble type (Ferrini & Galli, 1988). The corresponding 
efficiencies e M and en will be varied according to the pre- 
scription of Ferrini et al. (1994) while induced star forma- 
tion, due to the interaction of massive stars with the diffuse 
gas is considered as a local process, without any radial de- 
pendence, and its efficiency has the same value used in our 
previous disc models. Cloud formation and star formation by 
cloud-cloud collisions in bulge and core regions are enhanced 
in agreement with our previous analysis of disc zones. 

In the first application of the multiphase model to the 
bulge region (MF), this was considered as the natural con- 
tinuation of the disc, as moving radially from the Sun to the 
Galactic Center. Hence the parameters for star and cloud 
formation were calculated following their radial variation as 
in Ferrini et al. (1994). The first approach gave results in 
good agreement with observational data, showing that it 
does not appear necessary to invoke a merger or an accre- 
tion from extragalactic matter to explain the origin of the 
bulge. However, observational data show that the bulge is 
not the disc continuation toward the Galactic center, but it 
represents a distinct region, as is it apparent from IRAS and 
COBE results (Dwek et al. 1995) . In this paper we consider 
the galactic bulge as a distinct region respect to the disc, 
formed by the infall of gas from the halo over the bulge, 
without appreciable exchange of gas with the disc (Wyse & 
Gilmore 1992). We consider that the central galactic region 
is formed by three distinct zones and corresponding pop- 
ulations: the bulge, the halo over the bulge and the core. 
The bulge is modeled following Kent (1992) mass model as 
an oblate spheroid flattened on the galactic plane with semi 
major axis x = 2 kpc, semi minor axis z = 1.5 kpc, and with 
a mass Mb = 1-8 x 1O 1O M0. The core has the same shape 
of the bulge scaled to a semi major axis of 500 pc. The three 
zones interact by means of gas infall toward the center and 
through an outflow of gas driven by supernovae (SN) winds. 
The time scales for collapse from halo to bulge (th) and 
from bulge to core (tb) have been chosen to reproduce the 
correct final masses for the three zones. With this constraint, 
we have used the values th ~ 0.7 Gyr and tb ~ 10 Gyr. The 
collapse time we find in this way for the halo is shorter but 
not very different from that found in MF by extrapolating 
from the value in the solar neighborhood (th ~ 1 Gyr). As 
in MF, a SN-driven wind occurs when the thermal energy 
produced by combined SN explosions is larger than binding 
energy, calculated as Eb = GM(R)M ga3 /R, where M(R) is 
the total mass within radius R. 



© 1999 RAS, MNRAS 000, 



Bulges 3 



12.0 




T(Gyr) 



15.0 



Figure 1. SFR as a function of time for the three zones: bulge (solid line), halo (dotted line) and core (long dashed line). Line styles 
are the same for all figures. 



In comparison with MF we have updated a element 
nucleosynthesis from massive stars, in the mass interval 
8M < M < 1OOM . We use the computations of SN type 
II nucleosynthesis by Woosley & Weaver (1995) (WW95). 
WW95 calculate nucleosynthesis models for three different 
explosion energies and for several initial abundances. Initial 
abundances do not affect very much nucleosynthesis yields, 
while larger explosion energies increase them. We adopt 
their model with intermediate explosion energies (Model B) 
with initial solar abundances, to have a larger number of 
stellar masses. 

The model we are going to apply to the bulges of nine 
external galaxies is the same used for the Galactic bulge, 
but simplified since we do not include the core region. We 
assume a protogalaxy with a sphere of primordial gas whose 
total mass M to t is calculated using the rotation curve V(R): 
M(R) = 2.32 x 10 5 x Vrot(R) 2 x R (Lequeux, 1983). 



3 THE BULGE OF MILKY WAY GALAXY 
3.1 Star Formation history 

Figure 1 shows the Star Formation Rate (SFR) in the three 
zones. Star formation in the halo over the bulge is similar to 
the star formation in the halo over the solar region: it dom- 
inates at early times and is essentially switched off after 2 



Gyr. The time scale for star formation in the bulge is longer, 
approximately 5 Gyr, while the maximum of star formation 
is reached at t ~ 0.6 Gyr. According to our model there is a 
significant number of intermediate age (< 10 Gyr) stars in 
the galactic bulge, in agreement with observational data of 
Holtzmann et al. (1993). In the core the maximum appears 
at the same time, but the star formation rate reduces more 
slowly because of the residual infall of gas from the bulge. 
We can conclude that star formation starts in the external 
region and continues mainly toward the center after about 
0.5 Gyr; the bulge formation occurred after the formation of 
the halo, in partial overlap with the thick disc formation and 
almost completely before the thin disc produces efficiently 
stars. 



3.2 Element Abundances 

Figure 2a and 2b show respectively the time evolution of iron 
abundance, [Fe/H] and of oxygen abundances, [O/H], for the 
three zones. The enrichment goes on during the collapse, 
then the relation steepens as the radius decreases. The halo 
stars are very metal poor. The average metallicity [Fe/H] = 
— 1.49 is in good agreement with observational data by Zinn 
(1985) for halo globular clusters within a distance R — 3 

kpc ([Fe/H] = -1.47). 

The average metallicity [Fe/H] for the bulge stars agrees 
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Figure 2. Time evolution of (a) [Fe/H] and (b) [O/H] for the three zones. 



with the observational data by Mc William & Rich (1994 - 
hereinafter McWR94) for a sample of K giants in the BW, 
as it will be discussed in next section. 

In Table 1 we present the main features of star forma- 
tion in columns 2-4: SFR at 13 Gyr (SFR no ,„), the ratio of 
maximum to present SFR (S m /S n ), the present mass frac- 
tion in each zone. In columns 5-9 are shown the main fea- 
tures of gas enrichment: the time at which solar abundances 
are attained by iron and oxygen (pFe and to respectively), 
the present iron and oxygen abundances and the average 
metallicity (iron abundance). 

The abundance of a elements is of particular utility to 
understand the formation time scale of the galactic bulge. In 
Figures 3a-3d we show our model results for [O /Fe] , [Mg/Fe] , 
[Ca/Fe] and [Si/Fe] versus [Fe/H] for the three zones com- 
pared with the data by McWR94 for K giants in the BW. 
From the Figures 3 we can see that core stars present en- 
hanced abundances of a element compared with bulge ones. 
Observational data for O, Ca and Si are very well repro- 
duced. On the contrary, our model is not able to repro- 
duce observational abundance for Mg. The abundance ratio 
[Mg/Fe] found by McWR94 is enhanced by ~ 0.3 relative to 
solar over almost the entire [Fe/H] range. Only the two stars 
with the smallest [Fe/H] ratios are reproduced by our model. 
This underestimation is probably due to the SN II nucle- 
osynthesis model we have used. Thomas, Greggio & Ben- 
der (1998) compare SN II nucleosynthesis models by WW95 
with that of Thielemann, Nomoto & Hashimoto (1997), fo- 
cusing on the [Mg/Fe] ratios; they find that the last calcu- 
lation model leads to an enhancement of Mg abundance by 
20 per cent relative to WW95 models. To reproduce the Mg 
abundance, the SN II nucleosynthesis would be two times 
higher than predicted by WW95 model. We can conclude 
that no model available today can reproduce Mg abundance 
observed in the bulge K stars. 

Figure 4a shows the model result for Mg abundances 
compared with the data by Sadler, Rich & Terndrup (1996) 
(SRT96). The displayed data represent the average values of 



[Mg/Fe] computed dividing the stars sample into four bins. 
Model results are lower but still consistent with SRT96 ob- 
servational data. In Figure 4b we shows our model results 
for [(C+N)/Fe] ratios versus [Fe/H] compared with SRT96 
data, collected as in Figure 4a. The observational data seem 
well reproduced, even if the large uncertainty for abundance 
ratio (A[c+N/Fe\ = 0.4 dex) reduces the validity of this re- 
sult. 

3.3 Metallicity distributions 

The metallicity distribution of K giants is an important ob- 
servational constraint to understand the sequential time for- 
mation of the galactic bulge. Our model gives the time evo- 
lution of the total Star Formation Rate ^(t). To obtain the 
theoretical number of stars in the K giants phase from i/>(t) 
we used the isochrones from the FRANEC code and the 
Initial Mass Function adopted in this series of papers, de- 
rived by Ferrini, Palla & Penco (1990). All the isochrones 
correspond to an helium abundance Y — 0.27, according 
to Minniti (1996a), who finds Y = 0.28 ± 0.02 with a pho- 
tometric determination in the IR band of bulge stars. We 
used two different total metallicity: Z=0.01 for stars with 
age t < 2 Gyr and Z=0.02 for t > 2 Gyr. In Figure 5a- 
5c we show the theoretical K giants metallicity distribution 
for the three zones. All the distributions are normalized to 
unity. Moving from the halo to the core an increase of [Fe/H] 
and a decrease of standard deviation a[Fe/H] are evident. 
In Figure 6a we present the frequency histogram of Rich 
(1988) solution 1 [Fe/H] corrected by the regression relation 
of Mc William & Rich (1994) and the model result for the 
bulge region. Observations and theory agree reasonably in 
the range [Fe/H] > —1.25. For lower iron abundances the 
model predicts too many stars, as it happens with the G 
dwarf problem in the solar region. Pardi, Ferrini & Mat- 
teucci (1995) take into account the effect of the thick disc in 
the Galactic evolution; the presence of the this intermediate 
phase of accretion to the thin disc allows the solution of the 
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Table 1. SFR and metal enrichment 
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Figure 3. [o/Fe] versus [Fe/H] for the three zones, (a) oxygen, (ft) magnesium, (c) calcium (d) silicon. Data are from Mc William 
Rich (1994). 



G dwarf problem. We may easily imagine that the introduc- 
tion of a similar intermediate zone modifies the theoretical 
distribution on the low metallicity end. 

By comparing McWR94 observational data with theo- 
retical bulge K giants distribution of metallicity, we consider 
implicitly that every star in BW belongs to the bulge pop- 



ulation, without contaminations from halo and core pop- 
ulations. There are anyway indications of the existence of 
a contribution of halo star population in the BW: Rich 
(1990) finds a correlation between kinematics and metal- 
licity in a sample of K giants in the BW: giants with [Fe/H] 
< —0.3 have higher velocity dispersion than giants with 
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Figure 4. [Mg/Fe] (a) and (ft) [C+N/Fe] versus [Fe/H] for the three zones. Data are from SRT96. 



Table 2. K giants distributions: mean metallicity and standard 
deviations 
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- 1.60 
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0.45 
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- 0.60 
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[Fe/H]> —0.3. The velocity dispersion of the metal poor K 
giants is consistent with their belonging to the same r -3,5 
spheroid as globular cluster and RR Lyrae. Furthermore, in 
this paper we consider the existence of a third stellar pop- 
ulation, the core one, distinct from bulge population in its 
metallicity distribution. With this assumption a fraction of 
giants in the BW, inside 550 pc from the galactic center, 
belongs to core population. So we can conclude that in the 
BW there are not only present stars belonging to the bulge 
population, but there must exist an overlap of the three 
populations. 

In order to compare correctly theoretical results with 
observational data, we have to estimate the contributions 
from the different galactic components. To compute the con- 
tribution of the three stellar populations by varying the 
distance from Galactic Center we can use the correspond- 
ing mass distribution models. We have considered the mass 
model G2 of Dwek et al. (1995) for bulge star distribution; 
the halo star contribution has been modeled following an ax- 



isymmetric power law following Minniti et al. (1995). There 
are not available accurate core stars mass models; because 
of the stronger concentration, the contribution of core star 
population to the total star number can be relevant only in 
the BW region, up to 550 pc from the galactic center. So 
we have considered two combinations for the BW, differing 
in the core stars contribution. In combination 1, bulge stars 
represent 85%, halo stars 5% and core stars 10 % of the total 
stars. In combination 2 we have considered a larger contri- 
bution of core stars, choosing respectively 75%, 5% and 20%, 
for bulge, halo and core. Figure 6b shows McWR94 metallic- 
ity histogram compared with combination 1. Combination 1 
better reproduces observational data, even if for iron abun- 
dances lower than [Fe/H]=-1.25 our model still predicts too 
many stars. 

In Figures 7a and 7b, we show the metallicity distri- 
bution of pure bulge and combination 2 respectively com- 
pared with SRT96 data. SRT96 data are better reproduced 
by combination 2, in which core stars contribution is 20%. 
From Table 2 it is clear that the choice of combinations sat- 
isfies better the fitting of distributions. In Table 2 we insert 
also the properties of the combinations truncated to the low 
metallicity end, to simulate the presence of the thick disc. 

Figures 8a and 8b show Minniti et al. (1995) metallic- 
ity distribution of star samples in F588 bulge field at l,b — 
(8°, 7°) (8a) and in globular cluster M22 at I, b = (8°, 7°) 
(8b), corresponding to distances R = 1.5 — 1.7 kpc from 
Galactic Center. Observational distributions are compared 
with theoretical data of combination 3 (halo stars 30%, bulge 
stars 70%) which corresponds in our simple scheme to these 
bulge regions. The observational data are well reproduced 
by combination 3 in the metallicity range —2.0 <[Fe/H] 
< +0.4. It is still present an overestimation for star with 
[Fe/H] < —2. On the other hand the absence of theoreti- 
cal stars with [Fe/H]> +0.4 can be explained considering 
that the scatter in the metallicity observational determina- 
tions increases for the most metal rich giants. The fact that 
the errors increase with metallicity to A\F e /m — 0.4 for the 
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Figure 5. K giants metallicity distribution, (a) halo, (b) bulge, 
(c) core. 



more metal rich stars cast doubts on the reality of the values 
for the most metal rich stars. Indeed, it is not clear whether 
there are any stars in the samples with [Fe/H]> +0.5 even in 
the more central zone or core (Jablonka et al. 1996; Davidge, 
1998). 
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[Fe/H] 

Figure 6. K giants metallicity distribution by Mc William & Rich 
(1994) (dashed) compared with (a) bulge and (b) combination 1 
theoretical ones (solid). 



3.4 Metallicity gradient 

The determination of an abundance gradient in the central 
regions of the Milky Way is important for testing models 
of Galaxy formation. A metallicity gradient is predicted by 
dissipative collapse models. We have calculated the theoret- 
ical metallicity gradient comparing the average [Fe/H] of K 
giants distributions at different distances from the Galac- 
tic center. We use combination 1 for giants in the BW 
(R — 550 pc) and combination 3 to reproduce giants distri- 
bution at R — 1.5 kpc. We find the existence of a metallic- 
ity gradient d[Fe/H]/dr = —0.4 dex/kpc in the bulge region 
(0.5 < R < 1.5 kpc). We have calculated the metallicity gra- 
dient in the core region with the assumption that near the 
Galactic center, the core star population dominates. We find 
the existence of a metallicity gradient d[Fe/H]/dr= —0.8 
dex/kpc in the core region (R < 0.5 kpc), so that the neg- 
ative metallicity gradient is not limited to the bulge region, 
but it continues in the core region increasing in modulus. 
These values are in good agreement with the value found 
from spectroscopy by M95 and by the new IR data from 
FTK99 who find a value of -0.064 ± 0.012 dex/degree for 
the minor axis of our Bulge. 



© 1999 RAS, MNRAS 000, §-|l| 



8 M. Molld, F. Ferrini and G. Gozzi 




1.0 o.o 

[Fe/H] 



0.5 



0.0. 



JZ 



3.0 -2.5 -2.0 



-1.5 -1.0 -0.5 0.0 0.5 
[Fe/H] 



1.0 




1.0 0.0 
[Fe/H] 

Figure 7. K giants metallicity distribution by SRT96 (dashed 
line) compared with (a) bulge and (b) combination 2 theoretical 
ones. 
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Figure 8. K giants metallicity distributions by Minniti (1995) 
(dashed line) compared with combination 3. (a) globular cluster 
M22, (b) field F588. 



3.5 Gas and cloud phases 

We evaluate the time evolution of atomic and molecular 
gas for bulge and core zones respectively. At present time 
atomic and molecular gas are almost in the same quantity 
in the bulge region, while in the core region molecular gas 
exceeds atomic gas as from Table 3 where we compare the 
total abundances for atomic and molecular gas with obser- 
vational data from Sanders, Salomon & Scoville (1984) re- 
ferring to the zones with distances from the Galactic center 
R < 0.7 kpc and 0.7 < R < 1.5 kpc. The theoretical amount 
of molecular gas in the core region provides an underesti- 
mate to the observed quantity. On the other hand the total 
quantity of molecular gas estimated in both bulge and core 
regions together is smaller but in reasonable agreement with 
observational data. Our model takes into account infall of 
the atomic gas while molecular gas infall is not considered. 
The great concentration of molecular gas could only be re- 
produced if a similar infall parameter for the clouds in the 
bulge region is introduced. 

We must remind that a difference intervenes however 
from the nature of clouds: in the bulge region, clouds 
have larger densities and smaller dimensions (Blitz et al. 
1993) due to the gravitational potential and/or instabili- 



ties. The central half kpc region harbors a variety of phe- 
nomena unique in the whole galaxy. A number of dissipa- 
tive processes lead to concentration of gas into a "Central 
Molecular Zone" (CMZ) of about 200 pc radius, containing 
5 — 10 x 10 7 M Q of molecular gas. In the CMZ molecular 
medium is characterized by large densities (n > 10 4 cm -3 ) 
and a large filling factor (/ > 0.1). Such densities, usu- 
ally found only in the molecular clouds cores, are required 
to make clouds survive to tidal disruption. The clouds in 
CMZ also present high temperatures, in the interval 30 - 
200 K, typically ~ 70 K, large velocity dispersion (15-50 
km/s), and apparently high magnetic fields (for a complete 
review see Morris & Serabyn 1996). The bar potential and 
the strong tidal forces may enhance clouds velocity disper- 
sion, so that the dominant star formation process may be via 
relatively violent events, which cause strong compression of 
the already dense Galactic center clouds (Rich & Terndrup 
1997). 

psfile=fig9a.ps angle=270 hscale=40 vscale=40 
hoffset=-20 voffset=+250 
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Table 3. Gas Phases 



H 2 HI H 2 /HI 

(1O 7 M ) (1O 7 M ) 



bulge 2.3 2.50 0.9 

core 0.1 0.07 1.4 

R< 0.7 3 

0.7 <R< 1.5 ~ 1 

R< 1.5 3 - 5 0.10 50.0 



Table 4. Observational characteristics. 



Galaxy 


D 


T 


Arm 




V rot, max 


Name 


(Mpc) 




Class 


(kpc) 


(km/s) 


NGC 224 


0.7 


3 


12 


7.8 


250 


MWG 




4 




6.2 


220 


NGC 4303 


16.8 


4 


9 


4.6 


150 


NGC 4321 


16.8 


4 


12 


8.5 


270 


NGC 628 


7.2 


5 


9 


5.0 


220 


NGC 3198 


9.6 


5 




4.5 


160 


NGC 4535 


16.8 


5 


9 


9.1 


210 


NGC 598 


0.7 


6 


5 


2.6 


85 


NGC 6946 


5.7 


6 


9 


5.1 


180 


NGC 300 


1.6 


7 


5 


2.1 


80 



4.1 Theoretical Models 

Theoretical models have been calculated for bulges of nine 
spiral galaxies: NGC 224, NGC 628, NGC 3198, NGC 6946, 
NGC 598, NGC 300, NGC 4321, NGC 4303, and NGC 4535. 
These spiral galaxies have different Hubble type, implying 
that the corresponding bulges have different size and/or 
mass. Their contribution to the total light of the galaxy 
in the B-band is higher for early types and lower for late 
types. For later galaxies we will consider as the bulge the 
central region of the disc, even if it has not a clear bulge 
appearance. 

The model we apply is the same used for the Galactic 
bulge and described in section 2 but with a simplification: we 
do not consider the presence of the core component. Features 
for every galaxies used in our models are summarized in 
Table 4. For every galaxy, named in column (1), we have: 
the distance D in column (2), the Hubble type T in column 
(3), taken from Tully (1988) and Simien & de Vaucouleurs 
(1986); in column (4), the arm class, taken from Elmegreen 
& Elmegreen (1987) and Biviano et al. (1991); in column 
(5) we give the effective radius of the disc R e // in kpc and 
in column (6) the maximum rotational velocity. 

The total masses are obtained from the rotation curves 
referenced in MFD96 for six galaxies (NGC 224, NGC 628, 
NGC 3198, NGC 6946, NGC 598, and NGC 300) and in 
MHB99 for the three Virgo galaxies NGC 4321, NGC 4303 
and NGC 4535. Geometrical and mass data used as input 
in our models are shown in Table 5. Discs and bulges radii 
are in columns (2) and (3); the bulge and total masses are 
in columns (4) and (5). 

In Table 6, model input parameters are given: column 
(2) solar equivalent radius Ro, column (3) disc scale length 



Table 5. Input mass and geometry data. 







ixb 


Nig 




Name 


(kpcj 


(Kpc) 






NGC 224 


25 


4.0 


40 


4.35 


MWG 


20 


2.0 


18 


3.30 


NGC 4303 


11 


2.0 


8 


1.60 


NGC 4321 


16 


3.0 


.30 


4.00 


NGC 628 


16 


2.0 


10 


3.00 


NGC 3198 


12 


1.5 


7 


1.80 


NGC 4535 


16 


2.0 


9 


3.30 


NGC 598 


9 


0.5 


0.8 


0.50 


NGC 6946 


12 


1.0 


6 


2.25 


NGC 300 


7 


0.1 


0.1 


0.45 



Table 6. Model input parameters. 



Galaxy 


Ro 


R s 


T coll 






Name 


(kpc) 


(kpc) 


(Gyr) 






NGC 224 


10 


5.4 


3 


0.45 


0.50 


MWG 


8 


4.0 


4 


0.15 


0.08 


NGC 4303 


6 


3.0 


8 


0.20 


0.015 


NGC 4321 


9 


5.0 


1 


0.45 


0.10 


NGC 628 


8 


4.0 


1 


0.25 


0.01 


NGC 3198 


4 


2.7 


5 


0.20 


0.02 


NGC 4535 


6 


2.7 


5 


0.25 


0.02 


NGC 598 


3 


1.7 


8 


0.05 


0.005 


NGC 6946 


5 


4.5 


6 


0.18 


0.02 


NGC 300 


2 


1.7 


15 


0.07 


0.007 



R s , column (4) collapse time t(Ro); columns (5) and (6) 
e M and en, efficiencies respectively for cloud formation and 
cloud-cloud interaction, dependent on Hubble type. 

4.2 Star Formation Histories 

In this section we analyze the model results for the nine 
studied cases of bulges or central regions, as for late type 
galaxies. 

We start by analyzing the star formation related results, 
summarized in Table 7. As expected more massive and early 
type galaxies present a more intense initial episode of star 
formation, suggesting a mini-burst (see Elmegreen 1999) , as 
shown in figure 9. 

The time when the maximum of the star formation oc- 
curs -column (2) of Table 7- is approximately the same for 
all bulges, being between 0.4 and 0.9 Gyr, always before 1 
Gyr. Logically, the absolute value maximum of the star for- 
mation rate column (3), is larger for massive bulges 
than for those with lower total mass. Therefore, the stel- 
lar mass created in this early time step is larger for more 
massive bulges corresponding to earlier type galaxies. One 
important point is that differences between maximum and 
present values for the star formation rates - columns (3) and 
(4) of Table 6 - are much lower in panel b) than in panel a) 
of Figure 9. This means that the SFR is more continuous, 
with a more constant value in time for bulges of later type 
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Table 7. Bulge model results. 



( 1 o vi; 

VjrdlcLA. V 


+ 


* m 




A/T TT T 

Ml HI 


M TT 


M 


M „ /M 


N&1T16 




U vi 0/y r ; 


u vi 0/y r ^ 


fin 


fin 9m^1 


(a n 9 A/r^ ~i 
V lu lvi 0^ 




NGC 224 


0.40 


29.2 


0.38 


0.30 


0.30 


37.50 


0.05 


MWG 


0.60 


8.6 


0.22 


0.19 


0.20 


17.87 


0.09 


NGC 4303 


0.90 


2.4 


0.13 


0.12 


0.31 


6.88 


0.09 


NGC 4321 


0.45 


17.3 


0.30 


0.20 


0.39 


27.16 


0.08 


NGC 628 


0.75 


4.1 


0.14 


0.09 


0.44 


8.83 


0.07 


NGC 3198 


0.85 


1.9 


0.12 


0.08 


0.19 


5.98 


0.12 


NGC 4535 


0.65 


3.7 


0.13 


0.07 


0.30 


8.00 


0.09 


NGC 598 


0.75 


0.15 


0.02 


0.01 


0.03 


0.03 


0.50 


NGC 6946 


0.60 


1.5 


0.10 


0.03 


0.10 


4.62 


0.25 


NGC 300 


0.40 


0.05 


0.002 


0.0002 


0.0008 


0.024 


3.20 



galaxies. In any case, the enhancement of star formation for 
the bulge bursts are not very high in comparison with higher 
values reached in elliptical or starburst galaxies. 

The ratio between past and present values is correlated 
with the Hubble type or with the Arm Class, since the ratio 
^(m)/»l/(n) varies with continuity. Indeed the morphologi- 
cal classification of galaxies in Hubble types is related to the 
bulge - disc luminosities ratio; a correlation between Hubble 
type and SFR in the bulge is expected due to the existent 
correlation between total mass and Hubble type. We quan- 
tify this relationship, that must be reproduced by coherent 
models of galaxy evolution. 

From our model, early type galaxies have larger and 
more massive bulges and form a large number of stars in the 
first phase of the starburst that interests the bulge region. 
There also exist a correlation between the star formation 
rate in the bulge and the Arm Class of the disc; this means 
that the star formation of the bulge depends also on the 
dynamical features of the whole galaxy, as we underlined in 
previous papers where we analyzed the correlation of galac- 
tic model properties, their observational counterparts and 
the corresponding Arm Class classification scheme, viewed 
as directly related to large scale dynamical processes. 

4.3 Gas Masses 

The star formation processes consume the gas in the bulge; 
both phases are larger in absolute value of mass - columns 
(5) and (6) of Table 6 - in earlier type bulges, but this ef- 
fect is due to the higher values of total mass. Taking into 
account geometrical effects, we calculate surface densities 
for atomic and molecular gas; it results that bulges of later 
type maintain higher values for both atomic and molecular 
gas surface densities: earlier bulges consumed their gas more 
rapidly than the later ones. The total molecular gas quan- 
tity depends on two process: the consumption to form stars 
and the formation of molecular clouds from diffuse gas. In- 
deed, the ratio H2/HI has a peculiar behavior: it is lower for 
early type bulges because the molecular gas has been con- 
verted in stars, but it is also low for later bulges due to the 
slow process of molecular cloud formation. Only intermedi- 
ate type bulges, such those with T=5 if massive enough, have 
high values for the ratio H2/HI, which may results five times 
higher than in T=3 type. It is interesting to note that this 



trend is completely different for spiral discs, where observed 
H2 /HI quantity is larger for early type galaxies (Deveraux & 
Young, 1989), as also derived in MFD96. The same kind of 
behaviour occurs when the ratio H2/HI is represented ver- 
sus the Arm Class: only intermediate galaxies show higher 
proportions of molecular clouds in their bulges. 

The total stellar mass of bulge stars - given in column 
(7) of Table 6 - is decreasing with the Hubble type. For less 
massive bulges the collapse lasts a longer time, and a larger 
proportion of mass remains in the halo zone such it can be 
seen in column (8). Moreover the already collapsed mass 
needs a longer time to form stars because the efficiencies to 
form clouds and stars are lower. 

4.4 Gas Chemical Abundances 

We summarize our results concerning heavy element abun- 
dances in Table 8. A few aspects deserve a detailed discus- 
sion; in columns (3), (4), and (5) we show oxygen and iron 
abundances relative to hydrogen and magnesium over iron 
abundances at the time of SFR maximum. Since the bulk 
of stars created near that time, the average abundances of 
the corresponding bulge populations are reasonably well de- 
scribed by these values. Due to the differences in the star for- 
mation histories, we could imagine a priori that every bulge 
has a different metallicity. However, the total abundances at 
the SF maximum turn out very similar for all bulges, with 
values ~ Z0/IO. The oxygen abundances also are very sim- 
ilar, around —0.60 dex, while the iron abundance is lower, 
approximately around —0.75 dex (it is synthesized mostly 
by SN Type I, exploding with a delay with respect to the 
main episode of SF, and a dispersion exists depending on the 
total mass in the bulge). Therefore the oldest stellar gener- 
ation with the largest number of stars have approximately 
the same abundance for every bulge, independently of the 
Hubble type. 

For the present time, we have the results of columns 
(5) [0/H]„, (6) [Fe/H] n , (7) [Mg/Fe] n and (8)Z n . The total 
abundance almost reaches Zq in all cases with a weak cor- 
relation between Z n and the mass of the bulge, the Hubble 
type or the Arm Class. However there are differences be- 
tween elements. The oxygen abundance is about solar in all 
cases showing a saturation level for all bulges, but the iron 
abundance is higher for earlier Hubble types. 
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Figure 9. Star formation histories predicted for our bulge sample: a)for NGC 224, NGC 4321, MWG, NGC 618 and NGC 4535 bulges; 
b) NGC 4303, NGC 3198, NGC 6946, NGC 598 and NGC 300 bulges. 



The average abundances for stellar populations, calcu- 
lated taking into account all generations of stars are re- 
produced in Table 9 where column (2) is the average to- 
tal abundance Z, (3) average oxygen abundance [O/H], 
(4) average calcium abundance [Ca/H], (5) average iron 
abundance [Fe/H] and (6) average magnesium to iron ratio 



[Mg/Fe]. The stellar abundances are as expected interme- 
diate between maximum and present day gas abundances. 
The only possible comparison with data concerning MWG 
is certainly positive, since the mean stellar metallicity is 
[Fe/H] = —0.17; this value is in agreement with the most 
recently estimates from FTK99, who give a value of —0.2 
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Table 8. Chemical Evolution Results. 



Galaxy 


fO/Hl™ 

1 v -' / -M Til 


fFe/HU 


[Mg/Fel™ 

L J - v - l o/ Til 




\Fe/U]„ 


[Mg/Fel„ 




NGC 224 


-0.61 


-0.77 


0.33 


-0.01 


0.37 


-0.12 


0.029 


MWG 


-0.54 


-0.70 


0.33 


0.00 


0.31 


-0.05 


0.030 


NGC 4303 


-0.71 


-0.77 


0.23 


-0.02 


0.27 


-0.05 


0.025 


NGC 4321 


-0.59 


-0.76 


0.33 


0.00 


0.33 


-0.08 


0.028 


NGC 628 


-0.68 


-0.74 


0.23 


0.00 


0.30 


-0.06 


0.023 


NGC 3198 


-0.61 


-0.67 


0.23 


-0.04 


0.25 


-0.05 


0.024 


NGC 4535 


-0.67 


-0.71 


0.34 


-0.02 


0.29 


-0.06 


0.026 


NGC 598 


-0.78 


-0.77 


0.25 


-0.14 


0.11 


-0.03 


0.019 


NGC 6946 


-0.47 


-0.62 


0.28 


-0.03 


0.26 


-0.07 


0.025 


NGC 300 


-0.63 


-0.86 


0.40 


-0.12 


0.11 


-0.01 


0.019 



Table 9. Average Stellar Abundances. 



Galaxy 




z 


[O/H] 


[Ca/H] 


[Fe/H] 


[Mg/Fe] 


NGC 224 





.015 


-0.31 


-0.19 


-0.22 


0.11 


MWG 





.015 


-0.29 


-0.17 


-0.17 


0.08 


NGC 4303 





.014 


-0.36 


-0.24 


-0.24 


0.07 


NGC 4321 





.015 


-0.30 


-0.18 


-0.21 


0.10 


NGC 628 





.013 


-0.37 


-0.26 


-0.27 


0.09 


NGC 3198 





.014 


-0.31 


-0.19 


-0.17 


0.06 


NGC 4535 





.014 


-0.34 


-0.22 


-0.23 


0.08 


NGC 598 





.013 


-0.34 


-0.21 


-0.20 


0.06 


NGC 6946 





.016 


-0.24 


-0.11 


-0.09 


0.05 


NGC 300 





.016 


-0.23 


-0.10 


-0.07 


0.04 



dex and also with McWilliam & Rich (1994) (-0.25 dex) 
and the other author's data referenced in the Introduction. 
Metallicities predicted for our bulges show subsolar values 
for the stellar populations and hence seem to exclude any 
similarity to massive metal-rich elliptical galaxies. 

The different nucleosynthetic origin of a-elements and 
iron is evident in the behavior of [Mg/Fe] , higher than solar 
for older stellar populations. In Table 7 the value reached at 
t m is [Mg/Fe] ~ 0.2 dex, with a correlation with the Hub- 
ble type; [Mg/Fe] takes the highest values at t m , while at 
the present time it is below solar. At t m larger values are 
for earlier bulges, while at present larger values are for later 
bulges. The conspiracy of SF and nucleosynthesis from dif- 
ferent populations determine the average stellar abundance 
to follow a smooth variation with type. 



5 STELLAR POPULATIONS AND SPECTRAL 
INDICES 

Once evaluated the global evolution of the bulge region, in 
particular the history of star formation, and hence the se- 
quence of stellar generations and their corresponding ele- 
ment abundances, we may determine the spectral indices for 
the stellar populations, adopting the synthesis model results 
for the Single Stellar Populations (SSP). This method goes 
in the reverse direction of other authors (see e.g. Jablonka 
et al. 1996). 

We use the fitting functions given by Worthey et al. 



(1995) to assign the spectral indices Mg2, H/3 and < Fe > 
to each SSP, building blocks of the bulges, as functions of 
age and metallicity. We also assign a continuum flux, cal- 
culated with a new synthesis model developed by Molla et 
al. (1999), and finally we combine all generations in the pro- 
portion determined by the created stellar mass at every time 
step (or at every age). 

Results for the studied bulges are given in Table 10; in 
Figure 10 we show our predictions in the plane Mg2 — < 
Fe > with a set of spectral indices data observed in bulges 
by Jablonka et al. (1996), Idiart et al (1996a,b), Beauchamp 
& Hardy (1997), Vazdekis et al. (1997) and Samson et al. 
(1998), and the data from Gorgas et al. (1997) for low mass 
spheroidal galaxies. Our theoretical predictions are in agree- 
ment with the data observed in similar regions and with the- 
oretical model for SSP obtained by Worthey (1994) and by 
Borges et al. (1995) and Idiart & Freitas-Pacheco (1995). 
Spheroidal galaxies have spectral indices very similar to 
those of bulges. Galaxies with large values for both indices 
are out of the theoretical lines, probably in agreement with 
larger [Mg/Fe]. 

First we determined the star formation and chemical 
evolution, then the spectral properties, as a direct conse- 
quence of both. The advantages in affording the problem 
with this strategy are evident in comparison with the stan- 
dard analysis, based on the fitting procedure. In our study, 
from the first principles, late type galaxies form most stars 
when [Mg/Fe] has decreased in comparison with large bulges 
where stars are principally formed at early times, with high 
[Mg/Fe] values. Idiart et al. (1996b), to reproduce observa- 
tions, require [Mg/Fe] values sensibly higher than ours, with 
[Mg/Fe] reaching 0.6 dex. Our models have a mean stellar 
[Mg/Fe] between 0.11 and 0.05. The multiphase model gave 
a mean stellar bulge metallicity [Fe/H] ~ —0.2, in agreement 
with all recent estimations. If the statement of Peletier & 
Bacells (1996) that all bulges are similar in metallicities is 
true, our model is in agreement with them: our results of 
Table 9 give Z ~ 0.015 and [Fe/H] ~ -0.20 . Idiart et al. 
(1996a) for the Galactic bul ge requ ire [Fe/H] = -0.19, but 
Idiart et al. (1996b) obtain [Fe/H] = -0.02 (almost solar) 
for the bulges of external galaxies. This difference is striking 
and we cannot explain why those models result so different 
for the Galactic bulge and the others. In our case all bulges 
are similar in abundances. 

The spectral indices Mg2 and < Fe > show different de- 
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Figure 10. Spectral index Mg2 versus Fe52. 
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Table 10. Predicted Spectral Indices. 



Galaxy 


Mg 2 


Fe52 


H • 


name 


(mag) 


(A) 


(A) 


NGC 224 


0.19 


2.51 


2.32 


MWG 


0.19 


2.54 


2.41 


NGC 4303 


0.16 


2.36 


2.45 


NGC 4321 


0.18 


2.48 


2.05 


NGC 628 


0.17 


2.37 


2.38 


NGC 3198 


0.17 


2.43 


2.70 


NGC 4535 


0.17 


2.45 


2.33 


NGC 598 


0.15 


2.25 


3.11 


NGC 6946 


0.19 


2.54 


2.14 


NGC 300 


0.14 


2.13 


3.49 



pendence on the morphological type and/or the total mass 
of the bulge: Mg2 depends on T while < Fe > maintains 
constant. This results was already shown by Jablonka et al. 
(1996), when they represented their observed indices as a 
function of the luminosity M r . The luminosity has a corre- 
lation with the mass of the bulge, which in turn depends on 
the total mass of the galaxy and on the morphological type. 
Therefore, a clear dependence between the spectral index 
Mg2 and T is expected. The index < Fe > does not shown 
any trend with T, neither with M r . This different behavior 
is not explained by the different mean value [Mg/Fe] , neither 
by differences between iron abundances in the stellar popu- 
lations of different morphological types, but by the different 
mean ages of their populations. The earlier galaxies have old 
stellar populations with large magnesium abundances. Both 



effects increase the values of Mg2. Later type galaxies have 
young populations in larger proportions. These populations 
have low indices and the higher proportion in mass and lu- 
minosity in the integrated quantities decrease the values of 
indices Mg2 . For the index < Fe > the effect of the ages is 
the same, but iron abundances increase for later types, which 
produce more iron when stars are creating. Therefore, the 
integrated < Fe > is due to a mixing of old and metal-poor 
stellar populations and these young populations rich in iron 
in the later type galaxies, while in the earlier proceed only 
from old stellar populations. Both effect may produce sim- 
ilar results, and in consequence the final < Fe > does not 
depend on the morphological type or the total luminosity. 



6 DISCUSSION AND CONCLUSIONS 

The multiphase model, inserted in the dissipative collapse 
scenario following ELS picture, is able to reproduce the ob- 
served properties of bulges. The determination of the history 
of star formation as the principal ingredient of the evolution 
is useful in the understanding of the bulge variations in the 
Hubble sequence. In massive bulges corresponding to early 
type galaxies, the intensity of the star formation maximum 
is large enough to consume the whole gas, by preventing the 
subsequent star formation. The young stellar populations do 
not exist or they exist in very small proportions in compar- 
ison with the bulk of stars formed in earlier times. In later 
type galaxies, less massive, this maximum is not so large, 
and the star formation rate has almost a constant intensity. 
These results are in agreement with observations in the K- 
band from Seigar & James (1998a,b), whose results seem 
to indicate that the morphological types must be controlled 
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by the star formation ratio between past and present and 
not by the ratio bulge/disc, because it does not present any 
correlation with the morphological type T. 

Chemical results show bulges very similar in their abun- 
dances at the present time and for the time when the star 
formation maximum occurs, in agreement with estimations 
obtained from colors data (Peletier & Balcells, 1996; de- 
Jong, 1996a,b). Mean abundances have also been calculated 
by showing subsolar values and similar for all bulges in- 
dependently of the Hubble type, the arm class and/or the 
luminosity of their host galaxies. 

With a method reverse to the one adopted by Jablonka 
et al. (1996) and other authors, we compute the spectral 
indices Mg2 and < Fe > from our stellar generation re- 
sults, by using SSP synthesis models. Our results are in 
agreement with a set of observations for bulges in the plane 
Mg2— < Fe >. We may explain why the first of this indices 
depend on the Hubble type or the total luminosity of the 
bulge, while the index < Fe > does not show a similar trend. 
The evolution and the star formation histories of our bulge 
models are crucial for this respect: the mean age of the stel- 
lar populations in bulges are different, even if the formation 
of the bulge himself starts at the same moment. The star 
formation is more continuous in the central regions of the 
late type galaxies than in the early type ones, by producing 
more iron in the stars of the first ones. This mixing of the 
old metal poor populations and young and iron rich popu- 
lations in the late type galaxies produce the same result for 
< Fe > than is produced by the old metal poor population 
created in the early type galaxies in the first phases of star 
formation. The magnesium is produced instead very quickly 
in these last galaxies, by increasing the abundance of the 
stars. Therefore the index Mg2 show a clear trend with the 
morphological type because both ages and abundances have 
the same effect, by increasing the value of Mg2. 

The model shows that it is not needed the consider 
bars or radial flows to explain the features of bulges or their 
formation. Theoretical results for intermediate type bulges 
show the highest ratio H2/HI. These galaxies present the 
strongest observed bars, as from K band data; this suggests 
that bars may appear and disappear in the evolution of the 
bulge. They would contribute to produce a new starburst 
phase in bulges, but it may occurs several Gyr after the 
formation of the bulge. In this case, the gas flows into the 
central region, and younger stars will be created in the bulge. 
Their effects would be measured if the quantity of inflow gas 
is large enough or comparable to the mass of stars already 
created in the bulge. Otherwise, their influence will be neg- 
ligible in the indices Mg2 and < Fe > and we could only 
reveal them if we use other kind of data such the calcium 
triplet or H/3. 

We conclude presenting our answer to the questions ad- 
dressed in the introduction: 

• it is not needed to require the merging and/or accretion 
of external material to reproduce the main observed charac- 
teristics of bulges; the self-collapse of the protogalaxy seems 
to be sufficient to give the correct time scales, enrichment 
history and populations to reproduce observations not only 
for the Milky Way but also for external spiral galaxies. 

• The observed properties of bulge stellar populations are 
well reproduced in detail if a correct evaluation of the contri- 



bution from the various regions involved in the bulge forma- 
tion and evolution are considered, we mean the halo, bulge 
and core populations. 

• There is no evidence of an analogy with elliptical galax- 
ies: the apparent geometrical similarity does not correspond 
to any functional and physiological equivalence between 
these two families of astronomical objects. 

• The mean age of bulges varies according to the galac- 
tic type of the parent galaxy: there is a strong correlation 
between the bulge populations and the whole galactic evo- 
lution. 

In conclusion, not only it makes sense to figure bulge 
and disc evolution as two joint components of an intrinsically 
unitary structure, along the Hubble sequence, but it appears 
to be the only realistic way to obtain a good comprehension 
of bulge structure. 
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